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Abstract 

A partial wave analysis of pp — ?■ rjrjir^ data from the Crystal Barrel experiment is made 
in terms of s-channel resonances. The decay channels ao(980)r/, /o(1770)7r and /o(2105)7r 
provide evidence for two 1=1 J^^ = resonances. The first has mass M = 2360 it 25 
MeV and width T = 300t^[]° MeV, and the second M = 2070 ±35 MeV, F = 310t^|]° MeV. 
There is also evidence for a J^'^ = 2"+ state with M = 2005 ± 15 MeV and T = 200 ± 40 
MeV, decaying strongly to ao(980)7r. 



The present work is part of an analysis of pp annihilation with isospin 1 = 1 and charge 
conjugation number C = +1 in terms of s-channel resonances. A combined analysis of data 
from final states Stt", 7r°?7 and tt^t]' is reported separately [1]. Here we focus attention on 
annihilation to rjrjn^. In these data, signals are visible from final states aQ{980)r], /o(1500)7r, 
/o(1770)7r and /o(2105)7r. They carry no spin and may be analysed simply in terms of Legendre 
polynomials describing the production process. 

We find that the largest contributions to rjrjTr^ data arise from S-wave final states with = 
0~; those from 1+, 2~, 3+ and 4~ are somewhat smaller. Hence these data give a rather direct 
determination of contributions with quantum numbers 0~. In contrast, Sir^ data contain weak 0~ 
contributions from the final state /2(1270)7r with orbital angular momentum L = 2 in the decay 
process. There, the 0~ amplitude is hard to separate from larger L = 2 /2(1270)7r amphtudes 
with = 2' and 4". 

The rjriTT^ data treated here have been presented in two earlier publications [2,3]. They may 
be fitted with the channels 



pp ao(980)r/ (1) 

^ 02(1320)77 (2) 

-> 02(1660)77 (3) 

^ /o(1500)7r (4) 

^ /o(1770)7r (5) 

^ /o(2105)7r (6) 

^ /2(1980)7r (7) 

-> /2(1270)7r. (8) 



Fig. 1 shows Dalitz plots for data at all nine beam momenta. Figs. 2 and 3 show projections on 
to masses of nr] and 7777 combinations; histograms show the fit described below. There are clear 
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peaks in Fig. 3 due to /o(1500) — )■ 7777, and in Fig. 2 due to ao(980) and 02(1320) — ?■ nr]. The 
latter two are stronger than is immediately apparent from the figure, since the and 02 peaks 
originate from only one of the two TT^r] combinations, e.g. from iT^rji] the other combination iT^r]2 
produces a broad background when projected on to M(7r°?7i), as one sees from Dahtz plots. 

In Ref. [3], it was shown that a small but highly significant peak in the 7777 channel requires 
the presence of /o(1770) at beam momenta 900-1350 MeV/c. In Ref. [2], it was also shown that 
data at 1525-1940 MeV/c require a strong /o(2105) — ?■ rji] signal. That resonance has also been 
observed in pp r]r)[4] and was first identified in J/^' — )■ 7(47r) data [5], where /o(2105) — )■ aa, 
and a stands for the tttt S-wave amplitude. Ref. [2] also presented evidence for a broad 2+ signal 
in ?7?7 with M = 1980 ± 50 MeV, F = 500 ± 100 MeV. There is further evidence for this broad 
state in WA102 data on central production of An [6]. 

We now describe the partial wave formula, given by equn. (9) below, used to fit the data. 
This is the same as in the accompanying paper [1]. A full fit is made to production and decay 
of all channels (l)-(8) in terms of pp partial waves up to — Each partial wave has a 
distinctive dependence via relativistic tensor expressions Z{9,a, f3) on production angle 9 of each 
resonance in the centre of mass, and on its decay angles {a, (3) in the rest frame of the resonance. 
The energy dependence of each partial wave amplitude is expressed as the sum of up to two 
s-channel resonances plus a background. These resonances, described by simple Breit-Wigner 
functions of s with constant widths, are found to cluster in the mass ranges 1930-2070 MeV and 
2220-2360 MeV. The background, where required, is taken as the high energy tail of a resonance 
below the pp threshold. This parametrisation guarantees that partial wave amplitudes satisfy 
the necessary condition of analyticity, relating the energy dependence of magnitudes and phases. 
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Figure 2: Projections on to r]7r mass; histograms show the partial wave fit. Numbers in each 
panel indicate beam momenta in MeV/c. 



Blatt-Weisskopf centrifugal barrier factors -B^(p) and -Bi(g) are included in the partial wave 
amplitudes; Bii^p) incorporates the correct threshold dependence on momentum p in the pp 
channel for initial orbital angular momentum i, and Bilq) likewise describes the dependence on 
L and momentum q in the decay to the final state, eg. ao(980)?7. The radius of the barrier is set 
to 0.83 fm from the determination in Ref. [4]. The ao(980) amplitude for decay to particles 1 and 
2 is described by a Flatte formula F(si2) and other resonances are described by Breit-Wigner 
amplitudes. In summary, the full partial wave amplitude for channel (1), as an example, is given 
by 

/ = ^Be{p)B,iq) J^fZjt-M^^^M^^^ + ns2s)Z,s], (9) 

where Gj are complex coupling constants and the sum i is over s-channel resonances and back- 
ground. Fitted parameters are Gj, Mj and Fj. The factor 1/p accounts for the flux in the 
entrance pp channel, and ppp is the phase space for this channel, 2p/y/s. Near the pp threshold, 
the S-wave cross section is then proportional to l/p\/s, the well-known 1/v law. In Fig. 4 below, 
fitted cross sections will be shown multiplied by ps^^'^, so as to display the resonant behaviour 
free of kinematic factors. The phase space for the final state is accomodated in fitting the Dalitz 
plot. 

The parameters of each resonance are scanned in at least 7 steps over a range of typically 
±100 MeV. For all quantum numbers other than 0~ and 2~, masses and widths are determined 
less precisely by 777771° data than by 37r°, but are consistent with Ref. [1] within errors. The 
37r° data have statistics > 100, 000 events per momentum whereas present data have typically 
5000-9000 events at each momentum. 
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Figure 3: Projections on to rjr] mass; histograms show the partial wave fit. Numbers in each 
panel indicate beam momenta in MeV/c. 
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Table 1: Columns 3 and 4 show masses and widths of resonances fitted to 777777° data. Columns 
5 and 6 show masses and widths from weighted averages with fits to 37r°, 77^77 and 77^77'. 
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Columns 3 and 4 of Table 1 summarise masses and widths fitted to present data. Errors cover 
systematic variations between decay channels; they also cover variations depending on whether 
small components are included in the fit or arc omitted. Table 2 gives a quantitative picture 
of the significance of each component in the fit. It shows changes in log likelihood when each 
channel is dropped and the remainder are re-optimised. Our definition of log likelihood is such 
that it changes by 0.5 for a one standard deviation change in one parameter. Hence a change in 
log likelihood of 20 is rather significant (~ 5a, bearing in mind the number of fitted parameters). 
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Table 2: Changes in log likelihood when each resonance is dropped from the fit, and remaining 
contributions are re-optimised. In addition 772(2250) 02(1660)77 produces a change of 107 and 
02(2255) ^ /2(1980)77 a change of 363. 

The black circles on Fig. 4(a) show the integrated cross section divided by a factor 3, for 
comparison with those of individual channels. These individual contributions are lower because 
strong constructive interference between channels contributes positively to the integrated cross 
section. There is a sizeable cross section for the final state 02(1320)77, see the dotted curve of Fig. 
4(a). However, it does not yield new physics. This channel comes largely from pp triplet states 
with = 1"*", 2"*" and 3"'". Resonances in these partial waves are better determinined by 377° 
data, where statistics are very high and there are strong decays to /2( 1270)77. The contribution 
to present data from — 4+ is negligible, probably because of the L — 3 centrifugal barriers 
for both production and decay. There is a small (7%) intensity from /2(1270)77*', but it has little 
effect on other fitted amplitudes. We find no significant contribution from /2(1525) or /2(1565) 
at any momentum. 

For — 0~, three resonances are included. One, 77(1800), is below the pp threshold and 
may simulate physics background. Two further resonances at 2070 and 2360 MeV are needed 
and are sufficient to provide a good fit to the data. Figs. 4(b) and (c) show 0~ contributions to 
the cross section. 

We discuss the highest 0~ state first. It receives large contributions from /o(2105)77, shown 
by the dashed curve on Fig. 4(c), and from 00(980)77. The /o(2105) appears only in S-wave 
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Figure 4: Values of X = cross sections multiplied by ps^/^ for individual channels, (a) black 
circles show the overall cross section scaled by 1/3 for comparison with summed contributions 
to /o(1500)7r, 02(1320)7; and ao(980)?7 from all J^. (b) 0~ contributions: the full curve is 
the summed intensity from /o(1500)7r; remaining curves refer to ao(980)?7; the dashed curve is 
the combined intensity from all resonances, the dotted curve the coherent sum of 7r(1800) and 
7r(2070); the chain curve refers to 7r(2360). (c) and (d): 0~ and 2~ intensities from individual 
channels, (e): summed l"*" and 3"'" signals and that from 2+ — > [/2(1980)7r]L=i. (f), (g) and (h): 
variations of log likelihood InL with mass for 7r(2070), 7r(2360) and the low mass 1^2 Oq (980)77 
signal. 
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production and is very secure. The /o(1500)7r channel makes only a small contribution with 
= 0~ from both resonances at 2360 and 2070 McV, shown by the full curve of Fig. 4(b). 
Table 2 shows that (980)77 makes decisive contributions to both resonances at 2360 and 2070 
MeV. For the ao(980)?7 channel, there is large interference between the upper and lower 0~ 
resonances; the combined contribution is shown by the dashed curve of Fig. 4(b). Despite this 
interference, contributions from ao(980)r/ are very stable. This channel and /o(2105)7r agree 
closely on a mass M = 2355 ± 25 MeV for present data. The cross section for production of 
/o(2105)7r rises at high mass faster than phase space for that channel and requires production 
through 7r(2360). Fig. 4(g) shows the variation of log likelihood as the mass is varied. All fits 
with a variety of ingredients give masses in the range 2337 to 2377 MeV. This range is used to 
assess the systematic error, which is much larger than the statistical error. The fit to 37r'' data 
gives a higher but distinctly less accurate mass 2385 ± 45 MeV. The weighted mean of 2360 ± 25 
MeV fits both TjTjTT^ data and Stt^ well. The width is much less well determined: F — 300^50*^ 
MeV. 

At lower masses, there is a sizeable contribution from ao(980)?7 and also from /o(1770)7r, 
shown by the full curve of Fig. 4(c). In the previous analysis of Ref. [3], it was shown that 
/o(1770) production peaks in the momentum range 900 to 1200 MeV/c. Table 2 shows that 
/o(1770)7r makes a highly significant contribution of 324 to log likelihood. It is particularly 
useful, since it is again produced only through the S-wave and does not contribute to 7r(1800), 
which is too low in mass. The cross section for production of /o(1770)7r does not follow phase 
space for that channel, but requires production through the resonance at 2070 MeV; there is 
very little contribution from the upper 0~ state. 

From the present rirjir^ data, the optimum parameters of the lower resonance are M = 
2070 ± 35 MeV, F = 310l^|]° MeV, where errors are mostly systematic. Fig. 4(f) shows log 
likehhood against mass when the background is fitted by 7r(1800). However, the possibihty of 
alternative descriptions of the background introduces a systematic error of ±35 MeV in the mass 
determination. The 37r° data of Ref. [1] also require the presence of 7r(2070), but again suffer 
from interference with a background term. From those data, the optimum is at 2090 ± 65 MeV; 
the width is large and poorly determined, 285 ± 75 MeV. Both sets of data are well fitted with 
a mass of 2070 ± 35 MeV. 

The two 0~ states lie rather higher in mass than corresponding 7 = states [4] . These were 
observed at 2010^^^ and 2285 d= 20 MeV. In Ref. [4], all states were found to lie close to straight- 
line trajectories of against radial excitation number, with an average slope of 1.143 ± 0.013 
GeV^ per excitation. The spacing of the two 0~ states observed here is consistent with this 
empirical rule. One cannot compare accurately with 7r(1300) because of the large uncertainty 
in its mass. The 7r(1800) is consistent in mass with the required intermediate state, but there is 
evidence from its decay modes in favour of interpretation as a hybrid [7]. The VES group has 
also reported evidence in up for a 0~ state at 1750 MeV; this is a favoured decay mode for qq 
states [7] and may be the second radial excitation. 

We turn now to quantum numbers — 2~ . In an earlier publication on the channel rjrjn^Tr^ 
[8], we reported evidence for a 7r2 with M = 1880 ± 20 MeV, F = 255 ± 45 MeV, decaying 
dominantly to 02(1320)?]. Here, we again find a strong 02(1320)?] S-wave contribution shown by 
the full curve of Fig. 4(d); the 2~ — )■ [a2^]L=o amplitude is distinctive because of its P2(cosq;) 
dependence on the decay angle a for 02 (1320) — )■ rjir. The fitted mass and width for this 
amplitude are closely consistent with those for 7r2(1880). The contribution to /o(1500)7r, shown 
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Figure 5: Angular distributions for production of ao(980) for events in the rjn mass range 960 to 
1020 MeV; histograms show the partial wave fit. Numbers in each panel indicate beam momenta 
in MeV/c. 



by the dashed curve in Fig. 4(d), is consistent at low masses with the small branching ratio 
observed for 712(1 880) to that channel in Ref. [8]. 

However, we observe here an additional aQ{980)r] contribution with — 2~, considerably 
larger than that allowed in rjriir^Tr^ data of Ref. [8]. It would be conspicuous there, but is 
absent. Its fitted magnitude in present data is shown by the dotted curve in Fig. 4(d). Table 
2 shows that it makes a very large improvement in log likelihood, namely 531. Both ao(980)77 
and /o(1500)7r signals for — 2~ are clearly recognisable from their production amplitude 
P2(cos^), which interferes distinctively with the large 0~ amplitudes; this interference is clearly 
visible at the lowest momenta in the angular distribution against production angle 9 in Fig. 5. 

If the 2~ — )■ [ao^]L=2 channel is fitted freely, it optimises at M = 1990 ± 30 MeV with a 
width of 290 ± 60 MeV. The variation of log likelihood with mass is shown in Fig. 4(h). Errors 
quoted for mass and width allow for the possibility of moving an aorj contribution into 7r2(1880) 
consistent with the upper limit from Ref. [8]. If the mass and width are set to those of 7r2(1880), 
despite the evidence in Ref. [8] against decays of 7r2(1880) to ao{980)r), log likelihood is worse 
by 31. For two degrees of freedom, this is more than a 7a effect. 

There is further evidence for two distinct 7r2 states from the 37r° data. They require a strong 
/2(1270)7r amplitude with L — 2, consistent in mass and width with 7r2(1880). The /2(1270)7r 
amplitude with L = peaks at a higher mass 2020 ±17 McV with F = 165 ± 35 MeV. The 
width is considerably less than found here. However, both those data and rirjir^ are well fitted 
with a compromise mass M = 2005 ± 15 MeV and width 200 ±40 MeV, together with 712 (1880). 

We have earlier found similar evidence for two neighbouring I = = 2^ states, 772(1860) 
and 772(2030) [9]. The 772(1860) has been confirmed by WA102 [10]. Two neighbouring = 2" 
states with I — 1 and different decay modes are then a clear possibility. If the spacing in 
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mass squared follows the empirically observed 1.143 GeV^, the radial excitation of 712 (1670) is 
expected at 1985 ±20 MeV. Because two separate 712 candidates in a narrow mass range require 
one of them to be an intruder state (probably the predicted hybrid), and because both lie at 
the bottom of our available mass range, we acknowledge that confirmation is desirable. We 
present the evidence so that other experimental groups should be alert to the possibility of two 
separate states. Present data are consistent with 7r2(1880) accounting for the entire 02(1820)77 
and /o(1500)7r signals; however, we cannot rule out the possibility that these two channels are 
also fed in the present data partly by the second state at 2005 MeV. 

For other quantum numbers, the rjrj'n'^ data are consistent with resonances required by the 
analysis of 37r°. There is a small but highly significant S"*" peak at 2260 ± 50 MeV in the 
dashed curve in Fig. 4(e). Despite large errors for mass and width, this is additional evidence 
that a resonance exists in this mass range. Likewise, there is a distinct peak in the 2~ 
amplitude for production of a2(1660)7r, shown by the chain curve of Fig. 4(d); unfortunately, 
the determination of mass and width are poor, both here and in Stt" data. There is a distinct 
4~ signal in ao{980)ri; Table 2 shows that it produces a large improvement in log hkelihood of 
213. It optimises at 2255 ± 30 MeV in present data, very close to the value 2250 ± 15 MeV for 
37r° data, where there is a large 4~ signal. For = shown by the full curve in Fig. 4(e), 
there is a definite signal at the higher masses; this is evidence in favour of the state required 
in the analysis of 37r° data at 2270 MeV. At the bottom of the mass range, a second 1+ state 
gives a small improvement of 37 in log likelihood. This improvement is sufficient to require some 
additional low mass contribution, but is not sufficient to determine the mass and width of any 
possible resonance. Table 2 includes small contributions from 2"*" states at 1950, 2030 and 2175 
MeV; those states are required by 37r° data, but contributions to t/t^tt^ are too small to help 
determine masses and widths. 

A distinctive feature of the 777777° data from 1350 to 1940 MeV/c is that they require a strong 
contribution from a broad /2(1980) decaying to 7777, channel (7). It is clearly visible by eye in 
the angular distributions for decay to 7777, see Fig. 4 of Ref. [2]. The data required the curious 
property that the /2(1980) is produced almost purely with spin projection m' = ±1 along the 
beam direction in the final state. The present analysis confirms this result but clarifies the 
reason for the helicity dependence. 

The large contribution from this state is shown by the dotted curve on Fig. 4(e). It peaks 
at ~ 2265 MeV, and is well fitted by an 02 (2255) which appears strongly in the analysis of 
37r'', TT^f] and 7r*^77' data. We find that the broad /2(1980) is produced purely by this 02 (2255) 
decaying to /2(1980)7r with L = 1. It has a ratio of coupling constants g between ^^2 and 
amphtudes of r2 — g{^F2)/g(^P2) — —1.9 ± 0.4 for present data; this ratio agrees well with the 
value —2.13 ± 0.20 determined in Ref. [1]. With this value of r2, the final state will be almost 
purely m' = ±1 for the following reasons. Clebsch-Gordan coefficients for coupling of pp to ^^2 
and ^P2 are such that the initial state will be purely tt^ = if r2 = — y^7/2 ~ —1.9. The m = 
state decays purely to final states with m' = ±1, again because of Clebsch-Gordan coefficients. 
Thus, the curious property that the /2(1980) is produced almost purely with m' = ±1 seems to 
be a fortuitous consequence of the fact that the 02(2255) has r2 close to —1.9. 

Figs. 5 and 6 show production angular distributions for ao(980)7/ and /o(1500)7r, selecting 
events within one half-width of the resonance mass. All are fitted quite well. Production and 
decay angular distributions for 02(1320)7; are illustrated in our previous publication [3] and are 
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Figure 6: Angular distributions for production of /o(1500) in the overall centre of mass for 
events lying within ±60 MeV of the resonance; histograms show the partial wave fit. Numbers 
in each panel indicate beam momenta in MeV/c. 



also well fitted. 

In summary, a partial wave analysis of r/r/vr^ data gives masses and widths for s-channel 
resonances consistent with those found in the analysis of Stt" data, Ti^r] and tx^t]' [1]. For 
= 0~, the rjriir^ data give the best determination of the mass and width of the state at 
2360 MeV. An additional 0~ state at lower mass is also required; data from Stt" and rjrjn^ are 
both well fitted with an average mass of 2070 ± 35 MeV. There is definitely a strong — 2~ 
amplitude in the ao(980)r7 channel at low masses, much stronger than observed in Ref. [8] for 
712(1880). It suggests a second 2~ state at ~ 2005 MeV, but needs confirmation because both 
2~ states are near the bottom of the available mass range. 
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